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Abstract 

A Ring Imaging Cerenkov detector built for the BRAHMS experiment at the 
Brookhaven RHIC is described. This detector has a high index of refraction gas ra- 
diator. Cerenkov light is focused on a photo-multiplier based photon detector with 
a large spherical mirror. The combination of momentum and ring radius measure- 
ment provides particle identification from 2.5 GeV/c up to 35 GeV/c for pions and 
kaons and well above 40 GeV/c for protons during runs that had the radiator index 
of refraction set at n — 1 = 1700 x 10 -6 . 



1 Introduction 



The BRAHMS (Broad RAnge Hadron 
Magnetic Spectrometers) Collabo- 
ration is one of the four heavy ion 
experiments running currently at 
the Brookhaven RHIC. This experi- 
ment was designed to study particle 
production from colliding systems as 
varied as p+p, d+Au and Au+Au in 
a wide rapidity range that extends 
from to 4 units of rapidity ( max- 
imum beam rapidity: 5.4 ) These 
measurements are primarily based 
on the extraction of inclusive trans- 
verse momentum distributions of 
fully identified particles with moder- 
ate reach in transverse momentum. 



Particle identification at momenta 
lower than 7 GeV/c is done via 
the time-of-flight technique and the 
more challenging task of identifying 
high momentum particles (up and 
above 30 GeV/c) is tackled with a 
photo-multiplier based Ring Imaging 
Cerenkov detector (RICH) described 
in the present paper. This detector 
is the last element of a two section, 
movable spectrometer called For- 
ward Spectrometer (FS). This spec- 
trometer was designed to measure 
high momentum charged particles at 
angles with respect to the colliding 
beams that range from 30 to 2.3°. 
Tracking and momentum measure- 
ment is done with four dipole mag- 
nets, two TPCs in the first section of 
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the spectrometer where occupancy 
can be high, and three drift chamber 
modules in the back section where 
lower occupancy and higher resolu- 
tion provide precise measurement of 
momentum. 

The BRAHMS RICH shown in Fig. 
[1] has a nominal radiator length of 
150cm, and a 55cm x 70cm spher- 
ical mirror with a radius R = 3m. 
The mirror was manufactured 
slump formed glass substrate later 
polished to best commercial mirror 
surface quality. The polished surface 
was coated with Al and a protective 
layer of SiO. The mirror is rotated 
by 8° to keep the photon detector out 
of the spectrometer acceptance. The 
photon detector is a photo-multiplier 
(PMT) based system described in 
section 3. This array of PMTs is 
placed on the focal plane of the ro- 
tated mirror. More details about the 
BRAHMS experimental setup can 
be found in [I]. 



2 Design and construction 



Cerenkov light is a particular form 
of energy loss present whenever a 
charged particle moves through the 
bulk of a medium and its velocity is 
higher than the velocity of light in 
the medium c/n; this emitted light 
appears as a conical "retarded wake" 
of light that follows the particle. 

The energy loss associated to Cerenkov 
light emission is evaluated at dis- 
tances much greater than the typical 
atomic lengths of the material and 
it appears as radiated photons of 
frequency u whenever the velocity 
of the particle (3 attains values such 
that the dielectric constant of the 
material e(u) and the particle veloc- 
ity satisfy the relation (3 2 e{uj) > 1 
see REf . [2] . The energy radiated per 
unit length is known as the Frank- 
Tamm relation: 




Fig. 1. Photograph of the RICH in- 
stalled at the end of the BRAHMS FS 
spectrometer. The interaction region is 
located some ~ 20 m to the right of the 
picture. 
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The integrand between brackets is 
equal to sin 2 (6 l ) where 6 is the angle 
of the emitted photon with respect 
to the particle velocity. From these 
relations it is possible to extract an 
estimate of the amount of light gen- 
erated in a particular detector. The 
fact that Cerenkov light emission 
does not depend on a microscopic 
description of the medium implies 
that the emission is coherent, and 
it is another quantum effect visible 
in macroscopic systems. The index 
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of refraction of the radi ator medium 
written as: n = \J tiyj) can be mea- 
sured for the frequencies of interest. 
A set of such measurements can be 
found in [3] and [I]. 

The conical wave of Cerenkov light is 
transformed into a plane wave by re- 
flection on a spherical mirror. If the 
angles of incidence to the mirror are 
small, all photons with the same an- 
gle of incidence are reflected onto a 
point on the focal plane, (located at 
R/2 where R is the radius of the mir- 
ror) and as Cerenkov light is emitted 
uniformly in the azimuthal angle <fi 
(defined with respect to the particle 
velocity), the focus of Cerenkov light 
reflected by an spherical mirror is a 
ring centered around a point related 
to the angle of incidence of the par- 
ticle with respect to the axis of the 
mirror. Within the approximation of 
small angles, an spherical mirror fo- 
cusing Cerenkov light on a photon de- 
tector placed on its focal plane mea- 
sures the charged particles angle of 
incidence and their velocity. 

The number of photon-electrons orig- 
inating from a path of length L in the 
radiator volume is written as: 

N detected = NqL SU1 2 (9) 

where Nq (in units of (cm) -1 ) con- 
tains all the detector parameters: 

Nq = 370 £ m T quatzfdet J QpAIT^duj 

where e m is the reflectivity of the 
spherical mirror set as indepen- 
dent of frequency and equal to 90%. 



Qpmt is the quantum efficiency of 
the photomultiplier photo cathode. 
The PMTs used in the BRAHMS 
RICH have borosilicate windows 
that set their quantum efficiency 
Qpmt lower wave length cutoff at 
250 nm and bialkali photocathodes 
with a maximum in efficiency of 20%. 
The calculated integral of Qpmt for 
these PMTs is equal to 0.33 eV. A 
fraction of the light, estimated to be 
equal to 9.7%, is lost when the pho- 
tons traverse the 2.5 cm thick quartz 
window. The transmission of that 
window r quartz is thus set to 0.90. Fi- 
nally the photon detector has dead 
between PMTs that amount 
to an fiducial term f det equal to 0.7. 
The expected figure of merit of the 
BRAHMS RICH is then equal to: 
N = 69 cm- 1 . 

The focused light doesn't form a per- 
fect ring, chromatic and geometrical 
aberrations will distort it. The posi- 
tion resolution of the photon detector 
will also contribute to the spread of 
the photons away from a perfect cir- 
cle. The angular or radial resolution 
of a Ring Imaging Cerenkov detector 
is critical at high values of momen- 
tum where the bands of radius ver- 
sus momentum for different particles 
begin to merge. The width of these 
bands sets the particle discrimination 
power. As mentioned above, the res- 
olution of the detector will be writ- 
ten as the sum in quadrature of the 
chromatic aberrations and the single 
hit radius resolution contribution: 



Ar\ 2 /Ar\ 2 

T J \ T J chromatic 
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. T J geometric \ T / detector 



The contribution from chromatic 
aberrations (Ar /r) chromatic is present 
because the index of refraction in 
a dielectric depends on the fre- 
quency of the electromagnetic wave 



as n = Je(uj). Photons of different 
energy focus at different radii. For 
a particular detector sensitive to a 
finite band of wave lengths, the con- 
tribution of this distortion to the 
resolution of the detector can be 
simplified at high momentum values 
where f3 ~ 1 and the maximum angle 
of emission can be written as: [I] 




n + l)(n - 1) 



n 



the chromatic aberration can then be 
written as: 
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The radius of a single focused photon 
has a fractional spread produced by 
chromatic aberration: 



/Ar 



AO 



An 



An in these relations is set by the dy- 
namic range of the photon detector 
and can be evaluated from measured 
values. In the case of the BRAHMS 
RICH the photon detector is an array 
of multianode photo-multipliers with 
a FWHM quantum efficiency equal to 
2.7 eV (from 250nm to 517nm). The 
radiator is a mixture of two heavy 
fluorocarbon gases C$F 1Q and C 5 F U 
from 3M [5]. 

The index of refraction of the gas 
mixture was calculated using mea- 
surements performed in liquid phase, 
[5] and the Lorentz-Lorenz equation. 
Assuming an average index of refrac- 
tion within the dynamic range of the 
PMTs equal to n - 1 = 1900 x 10~ 6 , 
the chromatic aberration contribu- 
tion to the resolution in the measure- 
ment of a single photon is equal to: 

(—) = 0.0466 

V T / chrom 

The four pixels of each R7600 photo- 
multiplier have square shapes with 
1.1 cm on each side. For each of the 
measured photons, the error intro- 
duced in the ring radius measure- 
ment by the assumption that the 
photon intersects the detector plane 
at the center of the pixel is listed as: 

/Ar\ s 

\ T J detector 



V r J chrom n(n 2 - 1) 
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with s = 1.1 cm and maximum ring 
radius of 8.7 cm the contribution 
to the radius resolution is equal to 
0.036. 

To estimate the overall effect of the 
geometrical aberrations we have sim- 
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ulated the geometry of the detector 
and ray traced all Cerenkov pho- 
tons produced by charged pions with 
fixed momentum set to 19 GeV/c. 
Two cases are considered, in the first 
one, a set of pions move all along 
the axis of symmetry of the detector. 
The deviation from a perfect ring is 
shown in the a and b panels of Fig[2J 
Panel a shows the distribution of the 
distance between the photon inter- 
section with the detector plane and 
a nominal center of the ring. The az- 
imuthal dependence of the deviation 
from a perfect ring (for pions at 19 
GeV/c and n — 1 = 1960 ppm) is 
shown in panel b. These deviations 
are all the result of the 9 degree ro- 
tation of the spherical mirror and 
the different positions of the emitted 
Cerenkov photons along the track 
of the charged pion. The photons 
with the smallest deviations are the 
ones emitted close to the mirror. The 
overall effect of mirror rotation and 
photon emission position is small 
and transforms the rings into ellipses 
with the big axis along the x axis of 
the photon detector. 

A second set of charged pions with 
the same fixed momentum was used 
to extract an overall effect from ge- 
ometrical aberrations. This time the 
tracks belong to detected particles 
and give us the most realistic distri- 
butions in slopes and intersections 
with the radiator volume. The result 
of ray tracing the Cherenkov photons 
from this second set of pions is shown 
in panels c and d of Fig. [2J The an- 
gles of incidence into the RICH are 
small, the most distorted rings are 
produced by charged particles that 
fly close to the edges of the radiator 
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Fig. 2. Geometrical aberrations from 
ray tracing in the actual geometry of 
the detector. Panels a and b correspond 
to the first set of pion moving all along 
the axis of the detector. Panels c and 
d were obtained with measured tracks 
that have different angles of incidence 
as well as different entrance point loca- 
tions, this set shows the strongest ge- 
ometrical aberrations. For both sets of 
pions, the left panel shows the distribu- 
tion of radii, and the right ones display 
the relative deviation from a perfect cir- 
cle as a function of the azimuth angle. 

volume. We can thus quote a max- 
imum value for geometrical aberra- 
tions as (Ar/r) geometric = 0.025 even 
though the majority of the rings are 
contained in the bright band of panel 
d that corresponds to 0.7% Once a 
particle is identified it is posible to 
correct the geometrical aberrations 
in an iterative way. The present anal- 
ysis does not include this correction. 

A prototype of this detector was 
tested in one of the experimental 
halls at the BNL AGS. These stud- 
ies aimed at developing a photo- 
multiplier based photon detector 
in collaboration with Hamamatsu 
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Corp. From an earlier PMT version 
with 256 pixels that provided the 
first rings, the development work 
continued to reduce the amount of 
charge shared between neighboring 
pixels. The second PMT version had 
100 x 1cm 2 pixel arranged in a 10X10 
matrix. This device included an ad- 
ditional focusing stage between cath- 
ode and first dynode reducing the 
charge spread to 10% at the center of 
the next pixel [6j, but had poor pulse 
height resolution, and was deemed 
too difficult to manufacture. Finally 
we tested the R5900 four pixel PMT 
mounted in a compact metallic can. 
In order to achieve close packing af 
these tubes Hamamtsu produced the 
R7600 03 M4F that was selected to 
be used in this detector. Results ob- 
tained with the prototype can be 
found in Refs. [6] and [7] . 



3 The photon detector 

The photon detector is an array of 80 
four-pixel photomultipliers R7600- 
03-M4F p]. Each PMT has a single 
photocathode plane evaporated on a 
borosilicate window. Eleven dynodes 
amplify the electron extracted from 
the cathode and the total charge is 
collected in four independent anode 
plates. A focusing stage is placed 
between cathode and dynodes to 
reduced the charge shared between 
anodes. The complete system is en- 
cased in a flange-less metallic can 
held at the cathode voltage. 

Close packing of these PMTs was 
achieved by Hamamatsu engineers 
with the design of a matrix of 5 x 



4 biasing bases and sockets. Each 
matrix has two high voltage con- 
nections such that a single channel 
of a HV power supply can bias ten 
PMTs. The photo-multipliers placed 
in these matrices were selected to 
have similar gains and pulse height 
resolutions. The signal from each 
anode is routed out of each matrix 
through RG17 cables. The complete 
photon detector consists of four PMT 
matrices, and 70% of the total area 
is covered by active pixels. 

The output of these PMTs is fast and 
has sufficient pulse height resolution, 
(in average the single photo-electron 
peak appears as a shoulder of the 
distribution) but their gain is insuf- 
ficient to send the signal through 
some 100 meters of cable (RG58) to 
be integrated in an ADC with 50 fC 
coulomb resolution. This shortcom- 
ing of the PMT was compensated by 
the addition of one stage of amplifi- 
cation mounted right on the matrices 
of PMT bases. Fast charge ampli- 
fiers were AC coupled to the anode 
output through impedance match- 
ing 24f2 resistors to ground and a 
0.1/j,f capacitor. Fig. 3 shows one of 
the bias array with the five amplifier 
cards connected to it. 

The response of the R7600 PMTs 
was studied with a well collimated 
LED positioned in front of the photo- 
cathode with a two axis stepping 
motor system set to move in steps 
of 1 mm. The amount of light pro- 
duced with green LED is constant 
during the duration of the scan of 
one photo-multiplier and can thus 
be used to study the position depen- 
dence of the response of a particular 
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pixel. As the LED is 1 mm into the 
neighboring pixel it has 20% of the 
average pulse height in the center. 
Two mm further away that fraction 
drops to ~ 6% and is negligible at 
the center of the next pixel. 

The 320 signal cables are connected 
to a LeCroy 1885F Fastbus ADC af- 
ter suitable delay to allow for a trig- 
ger definition for each RHIC collision. 
The charge on each channel was inte- 
grated during a gate that lasted 120 
ns. 



Fig. 3. Photograph of 5 amplifier cards 
mounted on one of the 5X4 PMT base 
matrices. 



4 Filling procedure. 

Work performed on the prototype 
detector proved that the only filling 
procedure that guaranties the high- 
est concentration of the gas mixture 
in the radiator volume is the one 
that starts by evacuating the vol- 
ume. The drawback of this method 
is the need of a vessel designed to 
withstand a full atmosphere pres- 
sure differential. Once the radiator 
volume was evacuated, a boil off of 



C5F12 was introduced till the pres- 
sure reached 392 Torr (or 41.3% of 
the final mixture pressure of 1.25 at- 
mospheres). After that, C 4 F 10 is sent 
into the radiator volume till the final 
mixture pressure is reached. Once 
the filling is done, a small sample of 
the gas was used to measure the in- 
dex of refraction by counting fringes 
in a Young interferometer as one of 
the split beams of laser light goes 
through a volume that is filling with 
the gas sample, and the other passes 
through an equal length of vacuum. 

Fig.H]displays a portion of the fringes 
detected with a PIN diode to convey 
their good definition that makes pos- 
sible to make measurement of the in- 
dex of refraction with a resolution of 
one part in a million. 
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Fig. 4. Interference fringes detected 
with a PIN diode while a sample of the 
radiator gas is brought into one of the 
cavities of the interferometer while the 
other is kept at vacuum. 

The highest index of refraction 
achieved with this mixture at 1.25 at- 
mospheres was n — 1 = 2029 x 10 -6 . 
Later, the focus of the collaboration 
shifted to studies at higher pt values 
and the gas mixture and operating 
pressure were changed to reduce the 
index of refraction to lower values 
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(n - 1 = 1600 x 10~ 6 ). 



5 Data analysis 

Tracks reconstructed with the FS 
spectrometer tracking system are 
used to find the nominal center of 
rings. In case there are several tracks 
in the event, the loop over tracks is 
ordered such that it starts with the 
one with the highest momentum. 
Once a ring center is defined, the dis- 
tance from pixel center to ring center 
is calculated for all pixels that have 
charge above pedestal. The radius of 
a ring candidate is set as the average 
of those distances. The set of pixels 
is accepted as a ring if their number 
exceeds a minimum set by default 
to be equal to 4 and is at least 40% 
of an estimated number of photo- 
electrons with a figure of merit set 
low as Nq = 55. 

Pixels that were included in the his- 
togram are marked and will not be 
used in the search for another ring in 
events where there is more than one 
track in the RICH. A more detailed 
description of this algorithm can be 
found in [9] . 



6 Performance 

Fig. [6] is a composite of five field 
settings of the FS spectrometer, no 
effort was made to normalize the 
yields, the main purpose of this fig- 
ure is to show the remarkable mo- 
mentum range of this detector; it can 
identify electrons, muons and pions 
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Fig. 5. Two rings of Cerenkov light fo- 
cused on the photon detector. The big 
ring (red pixels online) has been pro- 
duced by a 20 GeV/c pion. The smaller 
ring (blue pixels online) was produced 
by a 17 GeV/c proton. 

at low momentum, kaons are well 
separated from pions (at the three 
standard deviation level) up to ~ 25 
GeV/c. The index of refraction of 
the radiator gas throughout the runs 
was equal to n — 1 = 1690 x 10~ 6 
and the spectrometer was placed 
at four degrees with respect to the 
beam line. The rings with "satu- 
rated" radii extracted from runs 
where the index of refraction was 
set to n — 1 = 1560 x 10~ 6 have an 
average of 38 ± 9 photo-electrons. 
The measured figure of merit of this 
detector is thus N = 81 ± 16cm' 1 . 

Particle identification with the RICH 
detector can be done with two inde- 
pendent and consistent methods, the 
first one is based on the difference 
between measured ring radii and the 
expected radius of a ring produced 
by a candidate particle. If such differ- 
ence falls within a set number of stan- 
dard deviations the measured parti- 
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a m 2 is the standard deviation of the 
mass square distribution. 
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Fig. 6. The radius of the Cerenkov rings 
produced by negative particles at 4 de- 
grees with respect to the beams in p+p 
collisions at y/s = 200 GeV/c, as a 
function of their momentum multiplied 
by the charge. Different magnetic field 
settings of the FS spectrometer are in- 
cluded in this figure. No effort is made 
to normalize the different data samples. 

cle identity is set to be the one of the 
candidate particle. This method re- 
quires the correct value of the index 
of refraction of the radiator gas ex- 
tracted previously from the data and 
stored in a run information database. 
This method includes tools described 
in section 6.2 to handle high momen- 
tum particles whenever the pion and 
kaon band start to overlap. The sec- 
ond particle identification method is 
based on the calculated mass using 
the momentum of the particle, the ra- 
dius of the Cerenkov ring and the in- 
dex of refraction of the radiator. The 
resolution of this calculated mass is 
momentum dependent and has con- 
tributions from the momentum reso- 
lution as well as the radius resolution. 
Fig. [7] shows the distribution of mass 
squared as a function of momentum. 
The particle identification is done in 
this particular case with a ±2<r m 2 cut 
indicated with dashed curves, where 
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Fig. 7. Mass-squared as a function of 
momentum. The dashed curves show 
the ±2<r m 2 cu t used by the mass based 
particle identification method. The red 
dashed curve shows the threshold for 
Cerenkov light emission as function of 
momentum. 

Monte Carlo simulations of the 
BRAHMS spectrometer, together 
with information extracted from 
data, show the high efficiency of 
this detector; Fig. M displays the 
efficiency as function of the ratio 
l/lthrsh where jthrsh is the value of 
the 7 factor at the particle threshold. 
The efficiency values shown in this 
figure were obtained with protons. 
The simulations show a ~ 4% ineffi- 
ciency due to interactions with ma- 
terial at the entrance of the RICH. 



6. 1 Relative radius resolution 



The particle discrimination power of 
this detector is set by the relative 
radius resolution at each momen- 
tum value. Fig. [9] shows the standard 
deviation obtained from gaussian 
fits to the distributions of the ratio 

{j'measured T calc) / ' 'measured-, where 
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Fig. 8. Efficiency near threshold calcu- 
lated using protons. The horizontal axis 
is a normalized value of the 7 of the 
particle 

r ca i c is the expected radius for rings 
produced by pions, kaons or protons. 

The horizontal axis displays the ve- 
locity of the particles shown as their 7 
factor. Pions are shown with red filled 
circles and above 7 ~ 40 have a con- 
stant relative radius resolution with 
a value as low as 1.2%, kaons shown 
with blue filled triangles and the anti- 
protons shown with filled star sym- 
bols have a worsening resolutions as 
their momentum diminishes. 



6.2 Particle identification at high 
momentum 



Fig. [6] shows that kaons can be sepa- 
rated from pions with simple cuts of 
the radius of the rings for momenta 
smaller than ~ 25GeV/c, and that 
protons are well separated beyond 
the momentum range of the figure. 
In order to extend the separation 
of kaons and pions to higher val- 
ues of momentum it is necessary to 
parametrize the uncorrected abun- 
dance of kaons and pions in small 
momentum bins in order to allocate 
probabilities or weights to events 
where the kaon band has started 
merging with the pion band. That 
parametrization was obtained at 
each full field spectrometer setting 
by fitting projections of narrow mo- 
mentum bands (500 MeV/c) onto 
the radius axis. The functional form 
used for these fits was the sum of two 
gaussians with equal widths, the free 
parameters of the fit were the nor- 
malizations, centroids and the com- 
mon width. Fig. [TU] shows the results 
of those fits. 



0.08 










0.07 




i 






0.06 










0.05 










i; 0.04 
— 




* 






< 0.03 










0.02 




















0.01 





















_n ni 















100 


200 



Fig. 9. A fit to the width of the pion 
band. 



Once the abundance of kaons with re- 
spect to pions is know for a particu- 
lar full field spectrometer it is possi- 
ble to assign probabilities to events 
that lie in the overlap of the kaon and 
pion bands. Fig. [11] shows the sepa- 
ration between kaons (blue hatched 
histogram centered around 8.1 cm) 
and the more abundant pions shown 
as the red histogram centered at 8.4 
cm) in a wide momentum bin (from 
30 to 44.6 GeV/c). Protons are also 
present in this figure but their dis- 
tribution is not gaussian because at 
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7 Conclusion 
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Fig. 10. The abundance of positive 
kaons with respect to pions in p+p col- 
lisions at high rapidity (y ~ 3) obtained 
from fits to projections of narrow mo- 
mentum bands onto the radius axis. 



these momenta, the rings radii are 
still changing fast as function of mo- 
mentum. 



The BRAHMS RICH described here 
has performed very well throughout 
the six years of data collection at 
RHIC. Its extended particle identifi- 
cation range has been instrumental 
in the study of particle production 
at high rapidity in several nuclear 
system that include p+p, d+Au and 
several heavy-ion systems. The high 
resolution of the radius measure- 
ment together with the simplicity 
of the photo-multiplier based pho- 
ton detector have made this detector 
the most important tool among the 
other particle identification counters 
in the BRAHMS setup. 
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that particle discrimination between 
kaons and pions is still possible even 
though their radius versus momentum 
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of probabilities is described in the text. 
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